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The decay mode B® — > Df K allows for one of the theoretically cleanest time dependent 
measurements of the CKM angle 7. This contribution reports the world best branching 
fraction of this decay relative to the Cabibbo-favoured mode B® — > D~ir + based on data 
sample of 0.37 fb _1 proton-proton collisions at ^fs — 7 TeV collected with the LHCb detector 

in 2011, resulting in BR(B° S -> Df K ± ) = (1.90 ± 0.12 s " 1 ' ± 0.13 ayat t° 1 1 l f3/fd ) x 10~ 4 . 



1 Motivation 

The least precise direct measured parameter of the unitary triangle is the angle 7. The high 
abundance of bb pairs, together with an excellent proper time resolution, an excellent particle 
identification and trigger capability to select hadronic final states, allows the LHCb experiment 
to determine this parameter through a time dependent analysis using the 5° DJA"* decay. 
Unlike the flavour-specific decay B® — > D~ir + , the Cabibbo-suppressed decay B® — > DJK^ 
proceeds through two different tree-level amplitudes of similar strength. 

These two decay amplitudes can have a large CP-violating interference via B® — B® mixing, 
allowing the determination of the CKM angle 7 with small theoretical uncertainties through the 
measurement of tagged and untagged time-dependent decay rates to both the D~K + and D+K~ 
final states H Although the B® — > DfK^ decay mode has been observed by the CDF ^ and 
BELLE'S collaborations, at present its branching fraction is known with an uncertainty around 
23%S1. Moreover, only the LHCb experiment has both the necessary decay time resolution and 
access to large enough signal yields to perform the time-dependent CP measurement. 



2 The LHCb experiment 

The LHCb detector^ is a single-arm forward spectrometer covering the pseudo-rapidity range 
2 < rj < 5, designed for studying particles containing b or c quarks. The detector includes a 
high-precision tracking system (silicon and straw tube technologies) and a dipole magnet with a 
bending power of about 4 Tm. The tracking system has a momentum resolution Ap/p that varies 
from 0.4% at 5 GeV to 0.6% at 100 GeV, an impact parameter resolution of 20/iin for tracks 
with high transverse momentum, and a decay time resolution of 50 fs. Charged hadrons are 
identified using two ring-imaging Cherenkov detectors. Calorimeter and muon systems provide 
the identification of photon, electron, hadron and muon candidates. 

The analysis is based on a sample of pp collisions corresponding to an integrated luminosity 
of 0.37 fb" 1 , collected at the LHC in 2011 at a centre-of-mass energy yfs = 7 TeV. 



3 Selection 



The channels considered as signal in this document are the decays B° — > D~tt + , B® — > 
D~tt + and the B® — > D^K^ 1 . These decays are all characterized by a similar topology and 
therefore the same trigger, stripping and offline selection are used to select them, minimizing 
the efficiency corrections. 

The LHCb trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage which applies full event reconstruction. 

The decays of B mesons can be distinguished from the background by using variables such 
as the pt and impact parameter x 2 of the B, D, and the final state particles with respect to the 
primary interaction. In addition, the vertex quality of the B and D candidates, the B lifetime, 
and the angle between the B momentum vector and the vector joining the B production and 
decay vertices are used in the selection. In order to remove charmless background a requirement 
in the flight distance x 2 of the Dj from the B® is applied^. 

Further suppression of combinatorial backgrounds is achieved using a gradient boosted de- 
cision tree technique^. The optimal working point is evaluated directly from a sub-sample 
of B® — > D~ir + events in data selected using particle identification and trigger requirements. 
The chosen figure of merit is the significance of the B® — > DJK^ signal, scaled according to 
the Cabibbo suppression relative to the B® — > D~ir + signal, with respect to the combinatorial 
background. Multiple candidates occur in about 2% of the events and in such cases a single 
candidate is selected at random. 

Particle identification (PID) criteria serve two purposes: separate the Cabibbo-favoured 
from the Cabibbo-suppressed modes (when applied to the bachelor particle) and suppress the 
misidentified backgrounds which have the same bachelor particle (when applied to the decay 
products of the Dj or D~). All PID criteria are based on the differences in log-likelihood (DLL) 
between the kaon, proton, or pion hypotheses. Their efficiencies are obtained from calibration 
samples of D* + -»• (D° -> K~ ~tt + )it + and A — > pir signals, which are themselves selected 
without any PID requirements. These samples are split according to the magnet polarity, binned 
in momentum and pt, and then reweighted to have the same momentum and p^ distributions 
as the signal decays under study. 

4 Mass fit 

The three signal decays are distinguished with particle identification requirements applied at the 
final stage of the analysis. The signal yields are obtained from extended maximum likelihood 
unbinned fits to the data. In order to achieve the highest sensitivity, the sample is fitted 
separately for the magnet up and down data. The signal line shapes are taken from simulated 
signal events. A mass constraint on the Dr s \ meson mass is used in order to improve the B mass 
resolution. 

The shape of the signal mass distribution is obtained fitting a double Crystal Ball function 
which consist of a common Gaussian with two exponential tails, one to describe the radiative 
tail present in the lower mass region and the other one describing the higher mass region where 
only the detector resolution is involved. 

A common signal shape describes properly both polarities so a simultaneous fit with a 
common mean and width of the double crystal ball function is used. The mean is free to float 
in all the fits, while the width is fixed in the i?!?-modes from the result obtained in data in the 
B° — > D~ir + fit corrected for the B° — B® differences observed in the simulation samples. The 
other parameters are fixed from simulation. 

Four sources of backgrounds are present: the combinatorial background, the charmless back- 
ground, the fully reconstructed (misidenfied) background and the partially reconstructed back- 



ground. The offline selection is optimized to reduce the combinatorial background contribution, 
and the remaining contamination is fitted with an exponential shape for the modes with a bach- 
elor pion, while it is taken flat for the B® — > DfK^ mode. The validity of this assumption is 
checked with wrong-side samples and accounted for in the systematic uncertainty associated to 
the fit model. The other two background categories have different components in the three fits 
and therefore are explained separately. In all the fits the partially reconstructed background 
shapes are obtained fitting a non-parametric function on samples of simulated events generated 
in the specific exclusive modes, corrected for the observed mass shifts, momentum spectra, and 
particle identification efficiencies observed in data when it is needed. The yields are left free 
when possible or a gaussian constraint is applied if an expected amount is computable. 

In the B° -> D~ -7r + mass fit the two relevant sources of partially reconstructed background 
are the B° — > D*~ir + and B° — > D~p + decays and their yields are left free to float in the fit. 

In the B® — > D~tt + mass fit the misidentified B° — > D~ir + shape is fixed from data using a 
reweighting procedure to account for misidentification momentum dependency^ . The number 
of expected events is computed from the yield obtained in the B° — > D~ir + fit and the PID effi- 
ciency obtained from calibration sample. Its yield is therefore constrained to this expected value 
with a 10% uncertainty. The B° — > D~ir + yield is calculated based on the B° — > D~ir + branch- 
ing fractional, the measured LHCb value of f s /fj^, an d the value of the B® — > D~ir + branching 
fractional. The shape used is the same of the signal B® — > D~tt + with the mean position fixed 
from the B° —> D~tt + fit. The partially reconstructed backgrounds relevant for this fit are the 
decays B® — > D*~ir + and the B® — > D~ p + . Due to the large correlation between these two com- 
ponents, a gaussian constraint is used for the fraction of these two backgrounds. The fraction is 
assumed to be the same as in the B° case, while the variation is assumed to correspond to 20% 
SU(3) breaking. 

In the B Q S -> Df K± mass fit there are numerous reflections which contribute to the mass 
distribution. The most important reflection is the misidentified B® — > D~ir + decay. Its shape 
is fixed from data using the reweighting procedure while the yield is left free to float. The 
same procedure is also applied on simulation sample to extract the shape of the B° — > D~K + 
misidentified background. The yield is constraint according to the expected B° — > D~K + 
yield corrected for the PID efficiency. In addition, there is potential cross-feed from partially 
reconstructed modes with a misidentified pion such as B® — > D~p + , as well as several small 
contributions from partially reconstructed backgrounds with similar mass shapes. The yields of 
these modes, whose branching fractions are known or can be estimated are constrained to values 
obtained based on criteria such as relative branching fractions and reconstruction efficiencies 
and PID probabilities^. The fit results are shown in Fig. Q] 



5 Systematic uncertainty 

Systematic uncertainties related to the fit are evaluated by generating large sets of simulated ex- 
periments using the nominal fit, and then fitting them with a model where certain parameters are 
varied. The sources of systematic uncertainty considered for the fit are signal widths, the slope 
of the combinatorial backgrounds, and constraints placed on specific backgrounds. The largest 
deviations are due to the signal widths and the fixed slope of the combinatorial background in the 
Bg — > D^K^ fit. The systematic uncertainty related to PID is evaluated using simulated signal 
and calibration samples. The observed signal yields are corrected by the difference observed in 
the (non-PID) selection efficiencies of different modes as measured from simulated events. A 
systematic uncertainty is assigned on the ratio to account for percent level differences between 
the data and the simulation. These are dominated by the simulation of the hardware trigger. A 
total systematic uncertainty of 3.9% for the ratio B® — > D^K^/B® — > D^tt + , of 3.4% for the 
ratio B® D--K+ B° D-tt + and of 4.6% for the ratio B° s Df K± B° D-tt+ is found. 




Figure 1: Mass distribution of the B° — > D~tv + candidates (top-left), B® — > DJ-k + candidates (top-right) and 
B" — > Df K (bottom). The stacked background shapes follow the same top-to-bottom order in the legend and 
in the plot. For illustration purposes the plot includes events from both magnet polarities. 



6 Results 

The sum of the B® — > D~K + and B® — > DfK~ branching fractions relative to B® — > Dji: + is 
obtained by correcting the raw signal yields for PID and selection efficiency differences and it 
leads to 

BR(B® Df K±) 

1 — rr = 0.0646 ± 0.0043 ± 0.0025 , (1) 

BR(B° S -^D7n+) V ; 

where the first uncertainty is statistical and the second is the total systematic uncertainty 

The relative yields of the three decays B° ->■ D~tt + , B® — > Djir + and B° s — > DJK 1 * 1 are 

used to extract the branching fraction of B® — > D~n + and B® — > DfK^ together with the 

recent f s / fd measurement from semileptonic decays^, leading to 

BR{B G S D-n + ) = (2.95 ± 0.05 ± 0.17±g;|l) x 10" 3 , (2) 
BR(B® — ^ DfK^) = (1.90 ± 0.12 ± 0.13i[J;Jl) x 10" 4 , (3) 

where the first uncertainty is statistical, the second is the experimental systematic plus the 
uncertainty arising from the B° — > D~ir + branching fraction, and the third is the uncertainty 
(statistical and systematic) from the semileptonic f s / fd measurement. Both measurements are 
significantly more precise than the existing world averages^. 
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